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CENTROID AND MOMENTS OF AN AREA USING A DIGITIZER 


by R. W. Patch 
Lewis Research Center 


INTRODUCTION 

The object of this program is to provide a general-purpose digitizer 
program for calculating area and first and second moments of any closed 
planar geometric figure such as a cross section of a beam or column or a 
spectral line shape that has been closed at the baseline. -For simple • 
figures such as circles , rectangles , or triangles , the area and moments 
are easily calculated by hand from simple formulas given in handbooks . 

For figures consisting of two or more simple figures superimposed, the 
area and moments may be calculated by hand with somewhat more difficulty. 
The principal application of this program, however, is to complicated or 
irregular figures such as spectral line shapes, beams, columns, rotating 
machine parts, statistical distributions, and ship hull sections, which 
are difficult or impossible to calculate by hand. 

In addition to the basic quantities of area and moments, certMn 
derived quantities are also calculated for convenience. 

Previous published digitizer programs have merely calculated area 
or circumference of closed geometric figures. 


PROBLEM TASK DESCRIPTION 

A closed planar geometric figure of arbitrary shape is shown graphi- 
cally in figure 1. Tlie size of the figure is such as to fit on the 
sensitive portion of the digitizer table (graphic tablet). The figure 
is to be digitized by going around it once clocki^ise with the cursor. 
Curved portions are to be digitized at closely spaced inteirvals (here- 
after referred to as "continuous digitizing" because the intervals are 
determined automatically). To improve accuracy, straight portions may 
be digitized by merely digitizing the two end points. 

The quantities to be calculated require several axes. Tlie digitizer 
axes X and y may be set anywhere on the sensitive portion of the digi- 
tizer table. Tlie closed figure has a centroid whose position is not known 
initially. Through the centroid are centre idal axes Xg and yg 
parallel to the x and y axes, respectively. Consequently, if the 
angular orientation of the x„ and yg axes relative to the figure is 
important, the figure should be rotated to the desired angular orienta- 
tion relative to the digitizer table before doing the digitizing for 
this program. This rotation may be facilitated by using the cursor and 
another small program to check the figure angular orientation (ref. 1). 

STAR category 09 
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The following quantities are to be calculated and displayed for 
each case. 

(1) Area of figure 

(2) First moment of x with respect to y axis 

(3) First moment of y with respect to x axis 

(4) Moment of inertia of area with respect to x axis 

(5) Moment of inertia of area with respect to y axis 

(6) Product of inertia of area with respect to x and y axes 

(7) X coordinate of centroid 

(8) y coordinate of centroid 

(9) Moment of inertia of area with respect to Xg axis 

(10) Moment of inertia of area with respect to y axis 

& 

(11) Product of inertia of area with respect to x„ and y„ axes 

O g 

(12) Polar moment of inertia of area around centroid 

(13) Radius of gyration about x axis 

(14) Radius of gyration about y axis 

(15) Variance in the x direction 

Most of these quantities are defined in reference 2, but they are also 
defined in, the next section or in the List of Symbols. 

METHOD OF SOLUTION 

Six basic integrals must be calculated. 

dA (1) 

X dA (2) 

y dA 



(3) 
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( 5 ) 


I 

xy 



xy dA 


( 6 ) 


(Symbols are defined in the List of Symbols.) From these integrals, all 
other required quantities may be found. To minimize storage requirements 
the six integrals are found by accumulating 12 sums as the closed figure 
is traversed with the cursor, assuming for curved portions that the fig- 
ure consists of a succession of straight lines between digitized points. 
The six integrals are calculated from the 12 sums at the end of the pro- 
gram. 


The area A may be expressed 


A = 



dy dx == 


s ^ 

upper 


dx 



(7) 


where "upper" and "lower" refer to the parts of the figure for which a 
perpendicular outward vector is in the +y and -y directions, respectively. 
Using the trapezoidal rule, equation (7) becomes 


A 


iZy ^^o + ^n^^^n - V 


( 8 ) 


where subscripts o and n refer to the first (old) and second (new) 
points, respectively, of an adjacent pair of digitized points, and the 
summation is over all pairs of adjacent points all the way around the 
closed figure. For Xn ^ ^o there is no contribution to equation (8). 


The first moment of x with respect to the y axis is 


= 




X dy dx = / xy dx - ^ xy dx 

upper lower 


( 9 ) 


or 


X 


- I Yj ‘■’‘n T Z/ '''■> ■ "o) 
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Again, there is no contribution if x = x„. 


The first moment of y with respect to the x axis is 





y dx dy ‘= / yx dy - / yx dy 

^ght left 


( 11 ) 


where "right" and "left" refer to the parts of the figure for which a 
perpendicular outward vector is in the +x and -x directions, respec- 
tively. Equation (11) becomes 


= I D 


- Vn>(yn 


-jYj ^o) 


( 12 ) 


There is no contribution to the last term when x = x . 

n 0 

The moment of Inertia of area with respect to the x axis is 




y^ dx dy = y^ x dy - y^ x 

^ right left 


dy (13) 


or 

I 


" i Xy yn^o + ^0) 

" ^ - ’^o^ ^^n yo^(yn + ^o) 

There is no contribution to the last term when x^^ = x^. 

The moment of inertia of area with respect to the y axis is 

2 


(14) 


I = 

y 


= ^y x^ y dx - ^ 

upper lower 


2 1 

X y dx 


(15) 


or 


" T ) , (y^x., " x„y^)[x^ + X X + x^) 
3 / / o n o-'n'\ n no 0 / 


? Z/ ‘’'n ' ’'o' "o) 
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There is no contribution if x, = x . 

n o 


The product of inertia of area with respect to the x and y axes 


is 


xy 


= n 1^ xy dy dx = ^ 

upper lower 


(17) 


or 


I = ^ 
xy 


‘S'..*!! - Vn> 


X^ - X 

n o 


+ 3 


^ (y„x - X y ) (y - y )(x^ + X X + x^^ 
1 \ •^o n o'^n n •'o'^x n no o/ 


n o 


1 Cyn- + 

+ 8 



(18) 


n o 


There is no contribution if x = x^ and, if this occurs, each rerm must 
be set equal to zero to avoid dividxng by zero. 

The other required quantities may be found from the six basic inte- 
grals. Tlxe X coordinate of the centroid is 




dA c 


X 

A 


(19) 


The y coordinate of the centroid is 




y dA Cy 
A " X 


( 20 ) 


The moment of inertia of area with respect to the Xg axis is (ref. 2) 


I = I - Ay^ 

gx X 


( 21 ) 


The moment of inertia of area with respect to the y axis is (ref. 2) 

O 


I = I - Ax 

gy y 


( 22 ) 
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The product of inertia of area with respect to the x and y axes is 
(ref. - 2) s S 


I = I - Axy 
gxy xy 

The polar moment of inertia of area around the centroid is (ref. 2) 

J = I ^ 

g gx gy 

The radius of gyration about the x axis is (ref. 2) 


(23) 


(24) 


(25) 


The radius of gyration about the y axis is (ref. 2) 


■V5 


The variance in the x direction is (ref. 3) 

iy 

®x A 


(26) 

(27) 


PROGRAM DESCRIPTION 

The program is written in BASIC for a Hewlett-Packard 9830A. calcu- 
lator and 9864A digitizer. 

The flow chart for MOMENT is given in figure 2, First an initial 
point XI, Y1 is digitized. Then a loop is entered which is traversed 
once for each new digitized point. During each traverse the sums T1 
and T2 are always added to , but the sums SO through S9 are only 
added to if X2 is not equal to X3. 

Closure occurs when the operator brings the cursor back sufficiently 
close to Che initial point. To check if this has occurred, two tests 
are made each time the loop is traversed: (1) has any point been outside 

the closure distance E from the initial point? (this test prevents 
closure before we have moved appreciably away from the initial point at 
the start), (2) is the last point inside the closure distance E from 
the initial point? E is arbitrarily set at 0.03 inch, and the dis- 
tances D from the initial point are approximated by 

^ ^il + ly^ - ^il < 28 ) 

when both conditions are met, closure is forced by making the new point 
the same as the initial point, and a last pass is made by the calculator 
through the upper part of the program loop. 
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When the last pass is complete, the loop is exited, and the digi- 
tizer gives two beeps to signal closure. The six Integrals are then cal- 
culated from the 12 sums, and the nine derived quantities are calculated. 
The six integrals and nine derived quantities are then displayed three at 
a time. 

The listing of MOMENT is 


10 

E = .03 


220 

Q = Y3 * Y3 + Y2 * Y2 

20 

PI = 0 


230 

N = Q + Y2 * Y3 

30 

F2 = 0 


240f 

T1 = T1 + B * G 

40 

ENTER (9,*) 

XI, Y1 

250 

T2 = T2 + G * N 

50 

X2 = XI 


260 

IP X2 = X3 THEN 420 

60 

Y2 = Y1 


270 

F = X3 - X2 

70 

S0 = 0 


280 

H = X3 + X2 

80 

SI = 0 


290 

I = Y3 - Y2 

90 

S2 == 0 


300 

M = X3 * X3 + X2 * X2 

100 

S3 = 0 


310 

L = M + X3 * X2 

110 

S4 = 0 


320 

S0 = S0 + B * F 

120 

S5 = 0 


33(7 

SI = SI + G * H 

130 

S6 = 0 


340 

S2 = S2 + I A L 

140 

S7 = 0 


350 

S3 = S3 + G A L 

150 

S8 = 0 


360 

S4==S4 + I*H*M 

160 

S9 = 0 


370 

S5 = S5 + F A N 

170 

T1 = 0 


380 

S6 = S6 + F A B A Q 

180 

T2 = 0 


390 

S7 = S7 + G A G A H/F 

190 

ENTER (9,*) 

X3,Y3 

400 

S8 = S8 + G A I A l/F 

200 

B = Y2 + Y3 


410 

S9=S9+I*I*H* M/F 

210 

G = Y2 * X3 

- X2 * Y3 

420 

IF F2 = 1 THEN 550 
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43Qf 

X2 

= X3 


44(7 

Y2 

= Y3 


45(7 

D = 

= ABSCX3 - 

■ XI) ■ 

46(^ 

IF 

D > E THEN 480 

47(7 

GOTO 450 


48(7 

FI 

= 1 


49(7 

IF 

D < E ANE 

> FI = 

5(7(7 

GOTO 190 


5117 

F2 

= 1 


52(7 

X3 

= XI 


53(7 

Y3 

= Yl 


54(7 

GOTO 200 


55(7 

WRITE (9,*) 


560 

WAIT 300 


570 

WRITE (9,*) 


580 

A = 

= S0/2 


590 

Cl 

= Sl/2 + 

S2/3 

600 

12 

= S3/3 + 

S4/4 

610 

C2 

= Tl/2 - 

S5/3 

620 

11 

= T2/3 - 

S6/4 

630 

13 

= S7/4 + 

S8/3 

640 

G1 

= Cl/A 


650 

G2 

= C2/A 


660 

14 

= 11 - A 

* G2 

670 

15 

= 12 - A 

* G1 


ABS(Y3 - Yl) 


1 THEN 51(^ 


S9/8 

G2 

G1 
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68(5 16 = 13 - A * G1 * G2 

690 J = I4 + 15 
700 K1 = SQR(I1/A) 

. 710 K2 = SQR(I2/A) 

720 V = 15/A 
730 FIXED 2 

74(/ DISP "XAVE"; Glj "VAR" j V; "A"; A 
750 STOP 

760 DISP "IX"; II; "lY"; 12; "IXY"; 13 
770 STOP 

780 DISP "IGX"; 14; "IGY" ; 15; "J"; J 
790 STOP 

800 DISP "IGXY"; l6; "CX"; Cl; "CY"; C2 
810 STOP 

820 DISP "KX"; Kl; "KY"; K2; "YAVE"; G2 
830 END 

The storage capacity required is 818 words. 

OPERATING INSTRUCTIONS 

Assuming the program is recorded on file 0 of a cassette, the fol- 
lowing is a check list for operating MOMENT: 

1. Turn on calculator and digitizer. 

2. Insert cassette. 

3. Line up axes of figure (if any) with axes of digitizer. 

4. Press SCATCH A EXECUTE. 

5. Press LOAD 0 EXECUTE. 

6. Set origin. 
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7. Press RUN EXECUTE. 


8. Place cursor at initial point of figure (anywhere on figure 
will do) , 

9. Press C on cursor to start continuous digitizing. 

10. Move cursor around curve clockwise back to initial point, 

11. Get t\jo beeps to indicate closure. 

12 •> Press C on cursor to turn off continuous digitizing. 

13. Read X, S^, A. 

14. Press CONT EXECUTE. 

15. Read 1,1,1 „ 

16. Press CONT EXECUTE. 

17. Read I , I , J , 

gx’ gy’ g' 


18. Press CONT EXECUTE. 

19. Read I , c , c , 

gxy» X’ y' 

20. Press CONT EXECUTE. 


21. Read k , k , y. 

x» y ’ 

22. Go back to step 6 or 7 for another case. 


If the figure contains a straight segment, approach the segment in the 
continuous C mode, press C to turn off the digitizer at the beginning of 
the segment, go to the other end of the segment by any path desired, 

press C to turn on the digitizer at the end of the segment, and proceed 

around the rest of the figure. The corners of a polygon may be digitized 

one at a time by pressing S once at each corner. 


A test case is given in figure 3, but should be laid out with draw- 
ing instruments. The correct results in appropriate powers of inches are 


XAVE 

1.00 

VAR 

0.67 

A 

8.00 

IX 

34.00 

lY 

13.33 

IXY 

14,67 

IGX 

7.11 

IGY 

5.33 

J 

12.44 

IGXY 

0.00 

CX 

8.00 

CY 

14.67 

KX 

2.06 

KY 

1.29 

YAVE 

1.83 





REPRODUCIBILITY 

OF . 





ORIGINAL PAGE IS POi... 
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However, do not expect to obtain accuracy to two decimal places on all 
quantities because the digitizer is only accurate to 0.01 inch and there 
is additional error in positioning the cursor. 
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APPENDIX - SYMBOLS 


BASIC 

Display 

Mathe- 

Description 



name 

name 

matical 

symbol 




A 

A 

A 

Area 



B 



y + y 

■'o ■'n 



Cl 

CX 

c 

X 

First moment of x with respect to y 

axis 

C2 

CY 

s 

First moment of y with respect to x 

axis 

D 


D 

Distance from initial point 



E 


E 

Closure distance from initial point 



F 



n o 



FI 



Flag to indicate if any point outside of closure 




distance 



F2 



Flag for last pass 



G 



y X - X y 
•^o n o''n 



G1 

XAVE 

X 

X coordinate of centroid 



G2 

YAVE 

J 

y coordinate of centroid 



H 



X + X 

n 0 



I 



^n ^o 
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IX 

I 

Moment of inertia of area with respect 

to 

X 




axis 



12 

lY 

I 

Moment of inertia of area with respect 

to 

y 



y 

axis 



13 

IXY 

I 

•VtT 

Product of inertia of area with respect 

to 

X 



xy 

and y axes 



14 

IGX 

H 

OQ 

Moment of inertia of area with respect 
axis 

to 


15 

IGY 

H 

00 

Moment of inertia of area with respect 
axis 

to 

"g 
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BASIC 

name 

Display 

name 

Mathe- 

matical 

symbol 


16 

IGXY 

I 

gxy 

Product of inertia of area with respect to 
and yg axes 

J 

J 


Polar iroment of inertia of area around the 
centroid 

K1 

BCX 

k 

X 

Radius of gyration about x axis 

K2 

KY 

k 

y 

Radius of gyration about y axis 

L 



2.2, 

X^ + X + X x 
n o no 

M 



2 u. 2 

n o 

N 



2 2 

y + y + y y 
■'n ■’^0 ■’^o-^n 

Q 



2.2 

so 




SI 




S2 



^(y„ - y„)(x^ + 

S3 



) " X y )(x^ + X X + x^') 

L / on o"^n ^ n no o/ 

S4 



Z/^yn - V^^n + V(^n-"^o) 

S5 



“ ^o>(yn + yn^o + ^o) 

S6 



Yy^r. - ^o^^yn + v(yn-^yo) 

S7 



V - ^oyn^’^^n + 

Lj 
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BASIC 

name 

58 

59 

T1 

T2 

V 

XI, Y1 
X2,Y2 
X3,Y3 


Display Mathe- 
name matical 
symbol 


+ Vo + \) 


n 


x_ 


VAR 


V ^oH^n ^o) 

/ ' X - X 

) ,(y X - X y ) (y + y ) 

/ / o n o-'n •'n ■'o' 

D (y X - X y )fy^ + y y + y^j 
■'o n o-^n‘'\-'n •'n-’^o •’^o/ 

2 

s Variance in the x direction 

^i’^i Coordinates of initial point 

XqjYq Coordinates of old point 

^n*^n Coordinates of new point 

x,y Digitizer coordinates or axes 

X jY Coordinates with respect to centroid or axes 
° ° through centroid 


T 
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FIGURE 1.- CLOSED PLANAR GEOAAETRIC FIGURE TO BE DIGITIZED, 
SHOWING DIGITIZER AXES AND CENTROIDAL AXES. THE FIGURE 
MAY BE ANY SHAPE. 



CALCULATE 6 INTEGRALS 
AND DERIVED QUANTITIES • 
^ I 

DISPUW Gl, V, A ) 

^ r ' 

DISPLAY II, 12, 13 J 

C DISPLAY 14, 15, J ) 

DISPLAY 16, Cl. C2 ') 

DISPLAY Kl, K2, G2 ) 

( ^~ END ') 


LAST PASS 



FORCE CLOSURE BY MAKING 
NEW POINT SAME AS 
INITIAL POINT 


FIGURE 2. -FLOWCHART FOR ’MOMENT'. 
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FIGURE 3.-TEST CASE (ALL DIMENSIONS IN INCHES). 



